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A system of two heavy fermions, leptons or quarks of the fourth generation, which are bound to- 
gether via the Higgs boson exchange is studied. The conventional Yukawa-type interaction produced 
by this exchange is accompanied by several important corrections. We derived the Hamiltonian, 
which describes the correction arising from the retardation (compare the Breit correction in QED); 
we also calculated the relativistic and radiative corrections. The Higgs-induced bound state appears 
for the fermion mass m > m cr ~ 500 GeV. When the long-range Coulomb interaction or the gluon 
exchange are included, the bound states exist for any mass, but the Higgs exchange drastically 
increases the binding energy of these states when m is approaching m cr . In the region m > m cr 
the gluon exchange gives a sizable correction to the Higgs induced binding energy. This correction 
greatly exceeds typical binding energies in the states produced via the gluon exchange only. The 
final results for the binding energies are presented as analytical functions of the Higgs and fermion 
masses. The possibility of detection of the considered bound states at LHC is discussed. 

PACS numbers: 12.39.Hg 14.80.Bn, 14.65.Ha, 
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I. INTRODUCTION 



Heavy fermions with masses 300-500 GeV and even 
heavier have been suggested recently to explain new re- 
sults on CP violation in B decays The bags with 
heavy fermions and bosons may play a crucial role in 
the bariogenesis @, Q . The nature of bound states made 
of such heavy particles is different from that for usual 
hadrons. The strength of the attraction between heavy 
particles due to the Higgs boson exchange increases with 
the particle mass. This may lead to formation of a new 
type of bound states. The Higgs-induced bags made of 
heavy fermions and bosons have been considered in nu- 
merous publications - see e.g. An alternative, 
but relative line of recent research addressed the problem 
of a possible dynamical electroweak symmetry breaking 
that may take place if the fourth generation is present 
28]. 

In several our previous works we concentrated mainly 
on "magic" bag consisting of 6 heavy quarks and 6 anti- 
quarks (altogether Y = 3x2x2 = 12 particles) occupy- 
ing Is shell. Its advantage is that it gives arguably the 
lowest limit on the fermion mass, which guarantees the 
binding. Besides, such multifermion bags may play a role 
in bariogenesis. However, these fermion bags are very 
difficult targets for possible observation at accelerators. 
Having this in mind in the present paper we consider the 
Higgs-induced bound states of only two fermions. The 
problem has recently been addressed in Ref. [29(, where 
the relativistic as well as radiative corrections induced by 
the fermion polarization were considered. To make our 
calculations accurate and reliable we consider a number 
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of important corrections including the relativistic and re- 
tardation corrections, as well as the radiative corrections. 
It should be emphasized that the retardation plays a very 
essential role in the problem, which makes it necessary 
to consider it alongside the relativistic correction. To de- 
scribe the joint contribution of the relativistic and retar- 
dation corrections we derive the nonrelativistic Hamilto- 
nian for the Higgs-induced potential, which includes all 
v 2 /c 2 relativistic corrections. This approach is similar 
in nature to the way the Breit correction is usually ac- 
counted for in QED. Our Hamiltonian-based treatment 
provides the convenient flexibility which allows it to be 
applied to a variety of similar problems in the future. 

To calculate the radiative corrections we employ the 
conventional methods developed previously to calculate 
the Lamb shift in atoms, see e.g. [3(| • Following the 
calculations of [I?} we express the radiative correction in 
terms of the renormalized Higgs vertex correction, which 
gives the fermion self-energy, and the renormalized Higgs 
polarization operator for the fermion loops. It should be 
stressed that the vertex correction derived for the first 
time in [27[ proves to be very large and has the sign 
opposite to the polarization correction, which makes it 
necessary to include it in the calculations. (Compare 
the Lamb shift problem, where the self-energy correc- 
tion dominates) . We solve the Schrodinger equation with 
an effective Hamiltonian which includes the nonrelativis- 
tic Yukawa Hamiltonian, the Brcit-typc relativistic terms 
and the radiative corrections, and find the ground state. 

Due to the short-range character of the interaction, 
the Higgs-induced bound state is formed only if the 
particle mass exceeds the limit, which is found to be 
m > m cr w (m ft /100GeV) 1/3 500GeV. When the long- 
range Coulomb interaction or the gluon exchange are in- 
cluded, the bound states exists for an arbitrary mass. 
We found that the Higgs-induced interaction produces 
an an interesting effect on these states even for m < m cr . 
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The Higgs induced interaction dramatically increases the 
binding energy and reduces the size of the bound state. 
For m — m cr the binding energy increases ah /a times 
and the size decreases ah /a times where ah is the con- 
stant of the Higgs exchange potential and a is the con- 
stant of the long-range potential. 



II. NON-RELATIVISTIC APPROXIMATION 

A. Yukawa potential 

It is convenient to start from the simple non-relativistic 
approximation. In the center of mass reference frame the 
interaction between two fermions due to the Higgs boson 
exchange may be described by the following Hamiltonian 

(h = c = l) 



H = + V Y (r) 
m 



(2.1) 



where Vy(r) is the Yukawa- type interaction induced by 
the Higgs exchange 



VV(r) 



- exp (~m h r) 



(2.2) 



Here m/, is the Higgs mass, and ah — g\l{^) is the effec- 
tive constant describing the strength of the interaction. 
Within the Standard model gh = m/v where v = 246 
GeV is the Higgs VEV. 

Numerical solution of the Schrodinger equation with 
the Hamiltonian (|2.1[) gives the bound state for 
ahm/2nih > 0.84. This corresponds to the particle mass 

m>m CT = 2v( 2.64 m h v 2 ) 1/3 w 2.76 (m h v 2 ) 1/3 . (2.3) 

A similar estimate was made previously in [28| . For rrih = 
100 GeV we obtain m cr « 503 GeV. This is larger than 
the critical mass m CIi bag = 320 GeV, which was found 
in for the bag of 12 fermions, but the distinction is 
not dramatic. To put it in perspective it is instructive to 
compare it with the situation when the bound state of the 
pair is produced by the heavy fermion interacting with 
some known, relatively light particle, when the critical 
mass proves to be much larger. For example, when the 
heavy fermion interacts with the t-quark the critical mass 
of the heavy particle is over 2 TeV. In the present work 
we do not consider so heavy particles. 



B. Adding Coulomb and gluon potentials 

The Higgs potential may be important even for m < 
m cr . Indeed, the bound state for m < m cr may be pro- 
duced by a relatively weak long-range interaction 



V h (r) 



(2.4) 



where a = a for the electrostatic interaction between 
charged leptons or 



— a, 
3 



(2.5) 



due to the gluon-exchange [31| . We restrict our attention 
here by the color singlet state of the pair. Generally 
speaking the pair of quarks can be in the octet state as 
well. However, it is known that the effective potential for 
the octet is repulsive, and what is worse, the octet pair 
of heavy quarks cannot exist by itself; its colored state 
needs the presence of additional, probably light quarks, 
which would compensate it. 

The constant a in Eqs. (l2.4[) . (|2.5I) is substantially 
smaller than a^,(m cr ) = 0.33. It is shown below that 
when the particle mass m is approaching m cr , the bind- 
ing energy increases ah/ a times. 

To estimate the energy we start from the variational 
approach, which gives accuracy about 10% (below we 
will perform accurate numerical calculations). Take the 
variational ground state wave function in the hydrogen- 
like form 



^)(r)=7r 1 / 2 q 3 / 2 exp (—qr) 



(2.6) 



Then from the Hamiltonian (|2.1[) one finds the following 
expectation value for the energy 



Eo(q) 



Aa h q 3 



m (m/j + 2q) 2 
The minimization of Eo(q) gives the equation on q 
dE {q) 2q 12a h q 2 Wa h q 3 



(2.7) 



m (m h + 2q) 2 (m h + 2q) 3 



= (2.8) 



Let us first find the critical condition for the bound state 
to appear. Combining equation (|2.8|) with E (q) = one 
finds 



<7cr = m h /2, 

(ah ™) cr /2 = ra h 



(2.9) 
(2.10) 



A similar estimate was made previously in |28| . The last 
condition allows one to find for the critical mass m cl - = 
2(irm h v 2 ) 1/3 w 2.93 (ra^ 2 ) 1/3 , which only slightly, by 
6% exceeds the corresponding value in (|2.3p extracted 
from numerical calculations. 

Let us add to the energy the contribution of the long- 
range potential Vj r = —a/r. Then we need to find the 
eigenvalue E of the Hamiltonian H — Hq + V^. Following 
the variational approximation we find then 



E(q) 
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dE(q) = 2q 
m 



Aa h q 3 



(m h + 2q) 2 
I2a h q 2 



aq 



(2.11) 



dq 



16ahq 3 
(m h + 2k) 2 (m h + 2k)'- 



- a = 



(2.12) 
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In the absence of the short-range potential, when ah = 0, 
these equations obviously reproduce the exact Coulomb- 
type wave function with q = ma/2 and energy E = Ec = 
— a 2 m/4. 

An interesting phenomenon takes place in the opposite 
case, when ah 3> a. Let us choose q = q cr = m,h/2, where 
q cl - is critical value, which was obtained in Eq. (|2.9p for the 
problem without the long range interaction. With this q 
the energy equals E = —aq = —arrih/2 = —a^am/ A, 
where Eq. (|2.10j) was used. Thus the Higgs exchange 
increases the binding energy of the system by a factor 
E j 'Ec — ah/ a. Correspondingly, the size of the system 
decreases by the same amount. We conclude that the 
Higgs contribution proves to be very important even well 
before the Higgs exchange alone would lead to the bind- 
ing. For leptons bound by the Coulomb attraction the 
Higgs exchange produces ~ 30 times enhancement of the 
binding near the critical mass. 



C. Confining potential 

The variational approach may also be used to estimate 
contribution of other interaction terms. In particular, in 
the region of large separation between the quarks the 
effect of the confinement could become essential. Note 
though that the region of large separations produces only 
small corrections to the energy. Therefore with the suf- 
ficient accuracy we can approximate an influence of the 
confinement using the simple approximation, by apply- 
ing the confining potential in its most transparent linear 
form V coni = or, where a = 0.2 GeV 2 [31]. We find then 
that the contribution to the energy equals SE colI = 2>a/q. 
For q > g cr = m.fc/2 this contribution is by three orders of 
magnitude smaller than the effect produced by the pure 
gluon exchange, which equals (4/3)a s q. 



Hence the effect produced by the Z-boson exchange is 
much smaller than the one originating from the Higgs bo- 
son exchange. One of the consequences is that a bound 
state of heavy neutrinos (from the hypothetical fourth 
generation) may be formed only if to > 500 GeV, and 
this bound state would be almost entirely due to the 
Higgs boson exchange. 



E. Higgs contribution for m < m CT /2 

If mass of the particle is significantly smaller than m cr 
the system is dominated by the long range interaction, 
i.e. q s» ma/2 -C rrc.fc/2, and the short-range Higgs ex- 
change is suppressed by the factor (ah/a)(ma/rrih) 2 ■ The 
Z-boson contribution has a similar suppression. 



F. Binding energy for m > m cr 

A relatively accurate result in the area of m > m cr 
may be obtained by using the expansion in powers of 
mh/2q. Indeed, near the minimum the energy E(q) is 
not sensitive to minor variations of q, and it is sufficient 
to find an approximate value of q. The error produced 
in this approximation of q may be smaller than the error 
of the variational approach itself. Keeping three leading 
terms in the mh/(2q) expansion of the Higgs contribution 
in Eqs. (|2.11|) , (|2.12[) we derive an approximate value 



q 



(a h + a) , 



(2.15) 



where ah is the effective value of the interaction strength 
suppressed by the short-range character of the Higgs ex- 
change, 



D. Z-boson contribution 

The weak interaction produces the short-range poten- 
tial via the Z-boson exchange, which with the sufficient 
for our purposes accuracy can be approximated as follows 



V z (r) 



a z 



exp {—mz r) 



(2.13) 



6th = olk ( 1 + 3 



m 2 (ah + a) 2 



(2.16) 



Substituting q from Eq. (|2.15[) into (|2.11l) we find the en- 
ergy, which for to > 0.75 TeV differs from the results 
of direct numerical solution of the Schrodinger equation 
only at the per cent level. 



Deriving this expression the low-energy, nonrelativistic 
scattering of the pair was considered in the unitary gauge. 
Using the variational approach one finds that f|2 . 13[) gives 
the following contribution to the energy 



Vz(q) 



Aa z q 3 



(m z + 2q) 2 



(2.14) 



For q > <7 cr = to/j/2 it is comparable to the contribu- 
tion due to the photon exchange (which can be approx- 
imated by the electrostatic interaction) since the short- 
range suppression factor (mz/2q + 1)~ 2 is close to 1. 



G. Binding energy in the limit m 3> 



In the limit to ^> the non-relativistic energy tends 



to 



E a fa -{a h + afm/A 



(2.17) 



However, the relativistic corrections ( oc afj become 
large in this limit since ah oc to 2 . For example, ah = 
m 2 /4irv 2 = 1 at m = 0.87TeV. In the following sections 
we will improve our results to account for this fact. 
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III. RELATIVISTIC CORRECTIONS 
A. Relativistic and retardation effects 

It is convenient to combine together the relativistic 
corrections for the propagation of fermions and the re- 
tardation correction for the propagation of the intermedi- 
ate Higgs boson. We rely on the conventional technique, 
similar to the one applied for derivation of the Breit cor- 
rections in QED. Consider first the scattering process for 
the two fermions interacting through the Higgs boson ex- 
change, see the diagram in Fig. Q] derive for this process 
the necessary relativistic and retardation corrections, and 
after that reformulate the result to make it applicable to 
the bound state problem. 





p. 


q, co | 




t 1 

p 2 ! 


p 2 



FIG. 1: Higgs boson exchange. The scattering process for 
two fermions is used in the text as a tool, which allows one 
to calculate the retardation and relativistic corrections. 

The analytical expression for the amplitude of the scat- 
tering process illustrated by Fig. [T]can be written as the 
matrix element 



M = (ui,U2 | V(q) | ui,u 2 ) 



(3.1) 



in which \ui) and (ttj| are spinors that describe the in- 
coming and outgoing fermions, while V represents the 
Higgs propagating in the intermediate state 



V(q) =gl 



fafa 



9l 



-q*- 

fafa 



(3.2) 



1 



Here fa, i = 1, 2 are the conventional Dirac matrices for 
the two fermions, which arise due to the scalar nature of 
the fermion-Higgs interaction (our notation for spinors 
presumes that (u = w). Considering the retardation as 
a correction we expanded the amplitude in powers of the 
transferred energy uj 2 . This energy and the transferred 
momentum q are related to the fermion energies and mo- 
menta ui = e[ - ex = £ 2 - 4) q = Pi Pi = Pa - p' 2 - 
As was mentioned, the term tx ui 2 in Eq. (|3.2[) repre- 
sents a correction that describes the retardation. Ne- 
glecting in this correction the relativistic effects we can 
use in this term the approximation fa w 1, i — 1, 2. We 
also can apply the nonrelativistic approximation relating 
u) to the fermion momenta 



Ul 



(Pi +q) 2 

2m 



P? 



P2 



(Pa - q) ; 

2m 



(3.3) 



Consequently we can rewrite (|3.2I) as follows 
fa fa (q 2 + 2p • q) 2 



,2\2 



(3.4) 



s q^ + mf i 4m 2 (q 2 -r 

Here we use the center of mass reference frame, p = 
Pi = — P2, which is natural for the bound state problem 
addresses below. 

The second term in the brackets in Eq. (|3.4[) is the con- 
tribution of the retardation to the relativistic correction. 
Let us simplify now the main, first term in the brack- 
ets in Eq. p.4p using the nonrelativistic approximation. 
With the help of the Foldy-Wouthuysen transformation 
we can relate the four-spinors u with the nonrelativistic 
two- spinors w 

a ■ p w 



(1 



i 

2m 



As a result we can use instead of the matrices fa 
following expressions 



(3.5) 

fa the 



fa 



1 



(p + q/2) 2 . (q X p) ■ ffj 



2m 2 



4m 2 



deriving after that 

fafa 

Here 



(p + q/2) 2 . (q x p) • s 
1 ^ I ■ 



2m 2 



S = - (<Ti 

2 1 



0-2) 



(3.6) 

(3.7) 
(3.8) 



is the operator of spin of the fermion pair. 
Combining Eqs.(pTl ]) .(|3"l) ]) we find 



V(q): 




2pq (qxp) 



2m 2 



4- 4m '(q-p) + 4 (q-p) 

4m 2 (q 2 + m 2 ) 



2m 2 



(3.9) 



In order to apply this result to the bound state problem 
we identify here p with an operator of the momentum 
p = iV and fulfill the Fourier transformation over the 
transferred momentum q. The coordinate r, which ap- 
pears as the result of this transformation, should be kept 
on the left side of the operator p. This machinery is very 
similar to the one, which is conventionally applied for the 
derivation of the Breit interaction in QED, see Section 83 
in 30] . As a result after straightforward calculations we 
find the following potential, which includes all relativistic 
and retardation corrections for two fermions interacting 
via the Higgs boson exchange 



Wrd = ST rcl + V Iei , 



P 



4m 3 ' 
2na h 



5{r) 



(3.10) 
(3.11) 
(3.12) 



an 



• exp (-to/, r) I -m h (1 — j + 



S + A 
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where 



3=1 



m h r \ d 2 



— + 1 - 2m h r + 



A 



2 'Or 2 
(l + m h r)l-s- 31(1 + 1) 
2r 2 



- — , (3.13) 

r or 

(3.14) 



It is presumed here that the state possesses the or- 
bital momentum 1 and spin s, and consequently the V- 
operator was transformed into expressions with d/dr. 
The first term ST rc \ in (|3 . 1 2|) describes the conventional 
rclativistic correction to the kinetic energy, while V ro \ 
originates from V(q). This term is written in a non- 
unitary form, which is similar to the way the Breit cor- 
rection is often presented in QED, though in both cases 
this fact produces no complications in the first order of 
the perturbation theory. However, if necessity arises 
the potential can be amended to make the unitarity 
transparent, V lci -)• V T ' el = (V lci + V? el )/2. Note that 
the term A vanishes for zero angular orbital momentum 
1 = 0. It is also interesting to note that the coeffi- 
cient ah/m 2 = (4irv 2 )~ 1 in front of 8(r) and exp (—nth r) 
terms in Eq. Q3.12p do not contain the mass m. 

Remember that deriving V(q) in (|3.9[) we treated the 
retardation correction, which stems from the second term 
in the brackets in (|3.2[) . by using the most simple non- 
relativistic approach and also rewrote the first term from 
the brackets in (|3.2[) using the Foldy-Wouthuysen trans- 
formation. However for some applications it may be use- 
ful to treat the problem using the relativistic four-spinors. 
The calculations presented in Appendix A show that with 
this purpose V(q) can be presented as follows 



-9l 



(7r(2p + q)) (-ft • (2p + q)) 



q 2 + m^ (cf+m 2 ,) 2 

(3.15) 

We verified that when the Foldy-Wouthuysen transfor- 
mation is fulfilled to replace the Dirac (3 and 7 matrices 
by their nonrelativistic approximations then Eq. (|3.15[) is 
reduced to Eq. (|3,9)l , Note though that an attempt to 
fulfill the Fourier transform over q directly in Eq. (|3.15[) 
should be treated with care since the matrix elements of 
7 matrices in p.ISp effectively depend on q. However, we 
will not dwell on this issue here since Eqs. (|3.9[) . (|3.12D and 
(13. 15)) constitute a reliable basis to treat the relativistic 
effects. 



B. Variational approach 

For the S state, in which 1 = 0, the last term in 
Eg. (13. 121) is absent, A = 0. Taking the variational wave 
function ip from (|2.6[) we find the expectation value for 
the relativistic correction from Eq. p.lOp 



(W re l) = (V|W r el|V) 



771" 



6a>hm 
2q + mh 



(3.16) 



Here the first and second terms in the brackets originate 
from the kinetic correction ST rc \ and the potential cor- 
rection VJ-ci respectively. The condition defining q and 
the corresponding nonrelativistic approximation for the 
energy E (q) were presented in Eqs. (|2.1ip . (I2.12[) . The 
total energy E can be found from E = Eo(q) + (W re i}- 

It is instructive to present this correction in relation to 
the kinetic energy T = q 2 jm 



rely 



T 



6a>hm 
2q + m h 



(3.17) 



C. Relativistic corrections to critical mass 

Taking the critical condition for the formation of the 
bound state from Eqs. (l2.9p . (|2.10[) we find that the right- 
hand side in (|3.17l) reads 



(W Tcl ) 



19 



m h 



4/3 



0.042 



(3.18) 



T 64tt 2 / 3 

Here the value = 100 GcV was adopted in the last 
identity. 

Obviously this small correction should result in a small 
variation of the critical mass. Let us find it, but first 
note that the considered relativistic correction represents 
only one of the different perturbations, which influence 
the critical mass. It is worth therefore to present an 
important simple formula for the shift of the critical mass 
in general case, for an arbitrary perturbation SV. Within 
the variational approach the necessary equation reads 



dm 



^(SV) _ (SV) 



3T 



(3.19) 



Here E (m) is the variational energy calculated without 
any corrections, while (SV) is the perturbation potential 
averaged over the variational wave function. To justify 
the first equality in Eq. (l3.19p one expands the condition 
for the critical mass Eo(m cr + 6m cr ) + (SV) = over 
Sm cr in the vicinity of the mass m cr , which satisfies the 
nonperturbed critical condition Eo(m cr ) = 0. To prove 
the last identity in Eq. (|3.19p one differentiates over m 
the energy Eq specified in Eq. (|2.7p . Remembering that 
ah = m 2 / (4irv 2 ) is a function of m, one finds 



dE 
dm 



m 2 m (mh + 2q) 2 



(3.20) 



where the last identity is derived from the critical con- 
dition En — 0. Thus one concludes that S-En = — 3T, 

u dm u 

which justifies the last identity in (|3.19p . 

Return now to the case of the relativistic correc- 
tion. Substituting the mentioned energy derivative into 
Eq. p,19p . specifying there that the perturbation is due 
to relativistic effects, 5V — > V IC \, we find the relative 
increase of the critical mass due to the relativistic and 
retardation corrections 



5m c 



(m 



reW 



3T 



0.014 



(3.21) 
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In the last identity here Eq. (|3.18p was used. In absolute 
units (5m cr ps 0.014m cr rs 7 GeV. 

We will see that the calculated in the next section ra- 
diative corrections reduce this mass shift by 30%, result- 
ing in the following total mass shift Sm cr ks 0.09to ps 5 
GeV. More accurate numerical results are mentioned be- 
low. 



D. Relativistic corrections for m S> mh 

The relativistic correction rapidly increases with mass 
77i. For large masses m 3> mh where q ps (6>h + a)m/2 
approaches q — cthm/2 — m 3 /(8nv 2 ) we obtain 



rcl 



T 



1.2a? 



3.4tj 



0.84 TeV 



(3.22) 



Thus, the naive perturbative treatment of the relativistic 
corrections based on W le \ from Eq. (|3.10[) definitely fails 
for m > 0.8 TeV. 

Moreover, the more accurate numerical treatment 
shows that this correction may be significant even at 
smaller masses. This perturbation has a positive sign, 
as can be picked up from the positive term oc 5(r) in 
Eq. p.l2[) . which gives significant contribution (~ 60%). 
The repulsion produced by this correction may be so 
strong that it destabilises the bound state. 

There is though the way around this difficulty. It will 
be discussed below, in Section IVI Al in detail. Speaking 
briefly, the idea is to allow the wave function the freedom 
to adjust itself in such a way that it is reduced at small 
separations r, where the relativistic correction is domi- 
nant. This diminishing of ip(0) reduces the relativistic 
correction and makes the perturbative approach sound, 
and the bound state stable. 

To prepare ourselves for this treatment of the relativis- 
tic correction we need to simplify Eq. (|3.12[) . With this 
purpose we, firstly, consider the most interesting for our 
discussion case of S-states, when the last term in this 
expression is absent, A = 0. Secondly, take into account 
that the Higgs mass is relatively small, mh <C m. As a 
result, in the region of small separation between the two 
fermions r < I/to the Higgs mass is negligible. Having 
this in mind we simplify the relativistic correction setting 
?77/j = in the term S, which gives p. 121) for 

d 2 1 d . 1 d ,„ no , 

S =TT + -V= A --7T- 3 - 23 
or z r or r or 

Here the last identity takes into account that we consider 
S'-states. Hence we find that for m/, = 1 = Eq. (|3.12j> 
gives 



rcl 



to V r or r 



to v ' r^ dr r 
To simplify notation the new parameter 



ro 



(3.24) 



(3.25) 



is introduced here (compare the classical radius of a 
fermion in QED, r c i as = a qc< i/m). To simplify ()3.24j) 
let us treat the sum of the first and second terms there 
as a perturbation. An inspiration comes from an identity 



1 d 



V>(r) 



dtp(r) 



dr = 0, 



27rr 2 dr T ' v ' J n r v ' dr 

(3.26) 

which implies that if one neglects these two terms from 
(|3.24j) and then contemplates the perturbation theory 
over them, one can be certain that the first order of 
this perturbation is absent. Neglecting the higher-order 
contribution of these terms (see discussion at the end of 
this subsection) we suppress these terms deriving from 
Ea. flXM 



ro 1 



(3.27) 



We can combine now this relativistic correction with the 
nonrelativistic Hamiltonian Hq from (|2.1I) and write the 
corresponding eigenvalue problem Etp = (Hq + V Te {)ip. 
Multiplying this differential equation by the factor r/(r+ 
ro) we present the result as follows 



■Eijj 





(El 

V TO 



Oth 



'•o 



(3.28) 



Clearly Eq. (|3.28j> can be considered as the eigenvalue 
problem of the Sturm-Liouville type. It has a very clear 
structure, but can be simplified further. Take into con- 
sideration that for not too large to, say m < 1 TeV, ro is 
expected to be smaller than the typical size of the bound 
state. Hence we can use the simplification rj(r + ro) « 1 
for the factor on the left-hand side of p.28[) . For large 
separation r ro the validity of this approximation is 
obvious. In the opposite case of small distances r < ro 
the term with E is much smaller then the potential en- 
ergy and can be ignored in the differential equation any- 
way. Thus the replacement r/(r + ro) ss 1 does not stir 
things up. 

We conclude that the wave function, which de- 
scribes the binding of two heavy fermions satisfies the 
Schrodinger- type eigenvalue problem E ip = i?o+reiV'i m 
which the Hamiltonian 



Hi 



Oth 



O+rel — 



'•() 



(3.29) 



accounts for the conventional nonrelativistic interaction 
and takes into account the relativistic and retardation 
corrections. 

A notable simplicity and clear physical nature of this 
result make it attractive. The only distinction of (13.29)) 
from the nonrelativistic Hamiltonian (|2.ip lies in the fac- 
tor ro, which effectively cuts the Yukawa-type potential 
off at small distances. As a result at small distances 
there arises an effective repulsion, which is in line with 
the repulsive nature of the relativistic correction men- 
tioned previously. 
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It is worth to summarize the set of approximations em- 
ployed for the derivation of Q3.29p . First, the Higgs mass 
was neglected in the relativistic term. Second, the ap- 
proximation r/(r + ro) « 1 on the left-hand side of (13.28)) 
was used. Both these approximations are justified by the 
parameters of the problem. Third, the sum of the first 
two terms in Eq. (|3.24p was neglected. We can see now 
the qualitative reason, which makes this approximation 
sound. The perturbation theory over these two terms 
starts from only the second order, as (|3.26p shows. Mean- 
while the repulsion at small separations caused by the 
relativistic correction is well accounted for in the Hamil- 
tonian iZo+rel- This repulsion reduces the wave function 
at small distances and hence suppresses the impact of 
the omitted terms in the problem. Consequently, their 
contribution is to be suppressed. This implies, in partic- 
ular that the sign of the correction to (|3.29[) is defined by 
the second order perturbation theory (over the omitted 
in (|3.24[) terms), and hence is definitely negative. 



IV. RADIATIVE CORRECTIONS 



To calculate the radiative corrections we use the con- 



ventional method which originally was developed to cal- 
culate the Lamb shift in atoms - see e.g. [3fJ. For the 
problem at hand the necessary corrections were calcu- 
lated in [221 . It was shown there that the most impor- 
tant contribution give two diagrams, shown in Fig. [5] 
One describes the Higgs polarization operator produced 
by the fermion loop, another the vertex correction to the 
Higgs-fermion interaction. 





(a) 



FIG. 2: Two most important radiative processes: (a) - the 
vacuum polarization produced by heavy fermions, which af- 
fects propagation of the Higgs boson, and (b) - the vertex 
correction for the Higgs-fermion interaction; solid and dashed 
lines - fermion and Higgs fields. 

Note that the diagram with the virtual fermion anni- 
hilation does not contribute to the energy shift for the 
non-relativistic fermions since the annihilation vertex of 
the fermion and anti-fermion into the Higgs vanishes for 
p 2 = 0. The reason is that the Higgs is a scalar. This 
is in contrast to a similar annihilation processes with the 
virtual annihilation of the positronium into photon, vec- 
tor particle, which gives a nonzero contribution in QED 
30}. We also neglect the diagram with crossed Higgs 
legs in the two-Higgs exchange since it is not enhanced 
by a large number of different fermions {y) in the loop 



(compare with the polarization operator Eq. (|4.ip 



A. Polarization operator 

In the approximation k <C m and m^ <C m the polar- 
ization operator represented by the diagram (a) in Fig. 
[5} equals [27| 



P{k 2 ) « -1 



9h ( fc2 " <f 
807r 2 m 2 



v (k 2 — ml) 2 

8^2 10^2 



(4.1) 



Here v is the number of heavy fermions with mass m 3> 
rrih (recall the factor N c = 3 for quarks). Importantly, 
the fermion mass m is canceled out in the final expression 
here since g^ = m/v. Hence the contributions of all heavy 
fermions are summed up. For one heavy lepton v = 1, 
for the whole 4-th generation v — 8, and v — 11 when the 
top quark is also counted. We neglect the contributions of 
W ± , Z° and Higgs in the loop since condition m ^> mt is 
not valid for them. The radiative correction to the Higgs 
potential Vh is given by 



P{k 2 )V h (k) 



g 2 (k 2 -m 2 )V h (k) 
807r 2 m 2 



(4.2) 



In the coordinate representation this gives a potential 
proportional to the (^-function (similar to the radiative 
correction to the energy in the positronium or hydrogen). 



B. Vertex 

Consider now diagram (b) in Fig|2j which describes 
correction T to the vertex of the fermion-Higgs interac- 
tion. The expansion of T(fc) in powers of k 2 /m 2 gives 



- T(k 2 
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1 



^U— = 1 + ^- (43) 
8tt 2 7 m 2 8tt 2 v 2 ' K 1 



where 7 is an expansion coefficient, which was calculated 
in [27]]. The found there analytical expression for 7 is 
lengthy, but a simple fitting proposed in this work, 



7 ~ 7fit 



I fl n m + mh 

3 V mft 



7m 



4m + 5?7jft 



(4.4) 



reproduces 7 with per cent accuracy. The correction to 
the Higgs-fermion vertex constant is Sgh = T(fc 2 ) — g/,. 
We find from this the correction to the Higgs exchange 
potential 



5V h ^k) = 2 -^V h (k) = ah7k2Vh{k) 



9h 



7T 777 



(4.5) 



C. Radiative potential and correction to energy 



Summing up the contributions of the polarization op- 
erator SV£ and the vertex SVF we find for them the 
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VWd(r) = ^(«h(7 ~ |j) *(r) + 7™. U h (r)) (4.6) 



following result, which is convenient to deal with using 
the coordinate representation 

Using the variational wave function ip( r ) from Eq. (j2.6p 
we obtain the radiative correction to the ground state 
energy 

4a 2 g 3 / v Airjml 



(V rad > 



V ' 20 



(m/, + 2q) : 



(4.7) 



For q = q CI = rrih/2 and the number of heavy fermions 
in the polarization loop v = 11 we obtain the ratio of the 
radiative correction to the kinetic energy 



(Vrad) 

T 



-0.015 



(4.8) 



We see that the radiative correction is about 3 times 
smaller than the relativistic correction Eq. p,18[) and is 
opposite in sign. For higher mass m there are significant 
cancellations between different terms in the radiative cor- 
rection Eq. (|4.7[) . As a result the radiative correction at 
to ~ 800 GeV is about 5 times smaller than the relativis- 
tic correction and has the opposite sign. 



V. VARIATIONAL ENERGY AT LARGE 
MASSES 

A. Relativistic kinetic energy 

We used previously the perturbation theory to account 
for the relativistic correction — p 4 /(4m 3 ) to the kinetic 
energy. This can be improved if we calculate explicitly 
the matrix element of the relativistic kinetic energy T re \ = 
2((m 2 -l-p 2 ) 1 / 2 — m) on the probing wave function (|2.6|) . 
which can be written as follows 

d 3 p 
(2*f 



(T rel ) 



( (TO2+p2) i/ 2 _ m y (p) ^. ( 5.i) 



where ip(p) — 871- 1 / 2 q 5 ^ 2 (q 2 +p 2 )~ 2 is the Fourier trans- 
form of the wave function (12 .6[) . Straightforward integra- 
tion yields 

<T rel ) 4 / x (3 - Ax 2 + 4x 4 ) (1 - 2x 2 )arccosxN 



7T V 



3(1 



r 2\2 



(l_ x 2)5/2 



(5-2) 

where x — q/m. The asymptotes of this expression 
are (T rcl }/m » x 2 (l - bx 2 /A + 128x 3 /15) when x -> 
and (T re \)/m s=s 16a;/(37r) — 2 for x —> 00. One veri- 
fies also that a very simple interpolating formula {T Te \) « 
m:r 2 /(l + 0.3x 2 ) reproduces the accurate expression for 
(T rc i) with errors below 12% provided q < 2m. Equa- 
tion (|5.2[) for the kinetic energy allows one to include 
the relativistic corrections into the variational procedure 
outlined below. This contrasts the naive perturbativc 
correction — p /4m 3 , which negative sign results in an 
unlimited increase of the momentum p during the mini- 
mization of the energy. 



B. Improved variational energy 

We can formulate now an improved version of the vari- 
ational energy, which incorporates the relativistic, retar- 
dation, long-range, and radiative corrections 

E(q) = (T rel ) + (V Y ) -aq+ (V ml ) + (V Iad ) . (5.3) 

Here all terms on the right-hand side are functions of q. 
The first term (T rc \) defined in (|5.2p represents the kinetic 
energy modified by the kinematic relativistic correction. 
The second one takes into account the Yukawa-type at- 
traction, 



(W) = -MK/r), 



Aa h q 3 



{m h + 2qf 



(5.4) 



being identical to the second term from Eq. (|2.7p . The 
term — aq comes from the long-range potential, gluon or 
photon exchange. The term (V^ i) incorporates the rel- 
ativistic corrections and the effect of retardation. Note 
that currently we include the kinematic relativistic cor- 
rection into (r rc i). In our previous arrangements the 
kinematic relativistic correction was taken into account 
in (ST Te {) via Eq. (|3.11|) . This correction shows itself 
as the first term in the brackets in Eq. (l3.16l) . (T re i) = 
— 5q 4 /(4m 3 ). We do not need to include this term here 
since it is covered by (T re \). The second term in the same 
brackets in (]3 . 16[) represents the correction to the poten- 
tial, which stems from the relativistic and retardation 
effects and contributes to Eq. (|5.3j) 



(V Tel ) = 



6a h q 4 



m 2 (2q + m/j) 



(5.5) 



The last term (V la .d) in (|5.3I) represents the radiative cor- 
rection calculated in Eq. (|4.7p . It is suppressed by a 
small factor ot^/n in comparison to the relativistic cor- 
rection. All five terms on the right-hand side of Eq. (|5.3[) 
are functions of the variational parameter q. Minimizing 
the total energy E(q) over q, one can find the variational 
estimate for the ground state of the two fermions. Nu- 
merical results are discussed below. 



VI. NUMERICAL DATA 

A. Perturbation theory 

To evaluate and study the behavior of the ground state 
energy for two heavy fermions bound together we employ 
the conventional approach based on the Schrodinger - 
type eigenvalue problem 



Etp(r) = Hip 



(6.1) 



We divide the total Hamlitonian into two parts, H = 
Hi + SV, and solve the Schrodinger equation (I6.ip accu- 
rately for H w Hi while applying the first order of the 
perturbation theory to the term SV . 
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In problems of this type the term SV often accounts 
for the relativistic and radiative corrections, which is the 
case in many atomic and nuclear problems. However, in 
the problem at hand we are facing a surprising devia- 
tion from this rule, the relativistic correction turns very 
large and repulsive. The repulsion is so strong that it is 
difficult to formulate the problem by choosing the naive 
Yukawa-type approximation, which ignores this repulsion 
at the initial step. If one tries to implement this naive 
approach and identifies Hi = H Q , where H a (|2.1I) in- 
cludes the pure Yukawa potential Vy(r) from (j2.2p . then 
the bound state found on the first step of such analy- 
ses would be destabilized, pushed into the continuum by 
the first order correction over the relativistic potential 

sv = w T ei (Euro . 

To overcome this difficulty let us have in mind that 
the relativistic correction is large predominantly in the 
region of small separation r between the fermions. The 
difficulty introduced by the naive Yukawa-type interac- 
tion is that it makes the wave function too large at the 
origin. As a result, the repulsion produced by the rela- 
tivistic corrections gains strength. 

The necessary remedy is clear. One needs to choose 
the initial approximation for Hi, which allows the wave 
function to be further extended in the region of large r, 
and correspondingly reduced at small r. Having this in 
mind we include into the Hamiltonian Hi the modified 
Yukawa interaction from f|3 . 29[) 

2 

Hi = Ho +I ei = ?h-e- m » r . (6.2) 

m r + ro 

Compared with the pure Yukawa interaction, its modi- 
fied version used here is not singular at the origin, and 
hence produces an effective repulsion at small r. Conse- 
quently it pushes the wave function out of the origin. An 
additional advantage of this approximation is that it has 
a clear physical meaning accounting for the relativistic 
correction, albeit in the simplified form. 

Eq. (|6.2|) neglects the long-range potential V\ r (r) in 
(|2.4j) . which stems from the possible gluon and photon 
exchange. The photon induced interaction should be con- 
sidered as a small perturbation. The neglect of the gluon 
exchange is obviously justified for leptons. For the quarks 
the problem is addressed in subsection IVI CI 

In the first order of the perturbation theory we sub- 
tract the approximate form for the interaction introduced 
in (|6.2|) . adding instead the necessary pure Yukawa- type 
potential and also adding the accurate relativistic poten- 
tial 

SV = V rcl res + V tad • (6.3) 

Kci, res = W TeX - a h (- ) e~ mh r (6.4) 

\r r + r J 

Here the term V^ e i, res can be called the residual rel- 
ativistic correction. It includes the proper relativis- 
tic correction W te \ from A3. 10)) and the additional term 
cx a exp (— wih t), which include the difference between 



the modified Yukawa interaction and the proper Yukawa 
potential. On top of it, the radiative correction Vr a d from 
Eq.(gl| is added into SV. 

Let us show now that Eqs. (|6.1[ )- (I6.4[) provide a sensible 
approach to the problem. To do this we calculated the 
binding energy e of the pair of heavy fermions. Fig. [3] 
presents this energy versus the mass of the fermion. The 
negative energy means that the bound state exists. The 
Higgs mass in this calculation was chosen rrih = 100 GeV, 
and the number of fermions contributing to the loop in 
the polarization operator was taken v = 11. 

Observe first of all that the way the perturbation the- 
ory is formulated Eqs. ([6.ip - (l6.4j) makes sense. Compare 
the dotted and thick dashed lines in FigJ3] The first one 
represents the eigenvalue Si derived from Eq. (l6.1[) when 
H — > Hi. The second shows e = e± + Se, where Se origi- 
nates from the residual relativistic correction K-ei.res from 
(|6.4p . Note a reasonable agreement between the two sets 
of data, their discrepancy is 29 — 38% for all the region 
shown, which makes the perturbation theory applicable. 
Hence, the accuracy of the found binding energy can be 
roughly estimated as ~ 15%. The found residual rela- 
tivistic correction is negative, which is due to the nega- 
tive sign of the second, additional term cx a exp {—rrih r) 
in Eq. (|6.4[) . This contrasts the positive sign produced by 
the proper relativistic correction V rc \ in this equation. 

Importantly, the bound state is definitely stable, it is 
not destabilized by the relativistic corrections. The cal- 
culations in Fig J3] demonstrate this fact in the first order 
of the perturbation theory. Moreover, if the second order 
of the perturbation theory over V rc i A - cs would be imple- 
mented (we do not fulfill these calculations, only con- 
template the expected outcome), then the negative con- 
tribution of the second-order correction would make the 
bound state only tighter. Thus, the way the perturba- 
tion theory is introduced makes the relativistic correction 
manageable. This is in contrast to the naive formulation 
of the problem, in which the pure Yukawa-type potential 
is taken as a starting point, while the relativistic correc- 
tion becomes a problem. 

The total ground state energy, which is shown in Fig. 
|3]by the thick solid line, includes also the radiative cor- 
rection calculated from V ra d in Eq. (|6.3[) . Observe that 
the latter is small and negative, compare the difference 
between the solid and dashed thick lines. 

Let us find the critical mass of fermions, which guar- 
antees an existence of the bound state. Fig. [H which 
shows the rescaled version of our data, allows us to state 
that m cr rs 515 GeV. Remember that Eq. (|2.3[) predicts 
7n cr Ps 503. It was mentioned that the pure Yukawa- 
based approach is not reliable, that the relativistic re- 
pulsion is strong. Nevertheless, we have to resign to the 
conclusion that the naive estimate of the critical mass 
based on the Yukawa- type approximation Eq. (|2.3|) is rea- 
sonably accurate. 

The discussed physical arguments and numerical re- 
sults are self-consistent, and hence look reliable. How- 
ever, we are to keep in mind that the problem considered 
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FIG. 3: Binding energy e of two quarks versus the fermion 
mass m; thick doted, dashed and solid lines are: the eigen- 
value of the Hamiltonian Hi (|6.2I) . the residual (see text) 
relativistic correction Se = (Vi-ci, res) (|6.4[) is added, the full 
correction Se = (8V), which includes the radiative correction 
is accounted for; thin lines: the variational energy e var from 
(|6.5(l with (solid) and without (dashed) the radiative correc- 
tion; m h = 100 GeV, v = 11. 




FIG. 4: The square-root of the binding energy (— e) 1//2 for 
two quarks versus the quark mass m, e is extracted from 
Eqs. (|6.1|) - (|6.4p and includes the relativistic and radiative cor- 
rections (same data is shown by thick solid line in Fig. [31); 
dots - calculations, solid line - interpolation. 



is unusual, the relativistic corrections are strong and re- 
pulsive. It is fortunate therefore that there is the way to 
verify our results by independent calculations discussed 
in the next subsection. 



B. Variational approach 

The wave function of the nonrelativistic problem has 
a very simple shape, possesses one maximum and shows 
no nodes. One should expect therefore that the sim- 
plest, one-parameter approximation for this function in 
Eq. (|2.6[) provides a steady basis for the variational ap- 
proximation. Following this path we considered the en- 
ergy E(q) found in Eq. (|5.3p as a function of the parame- 



ter q introduced in the wave function (12.61) . Minimizing 
this function one finds the binding energy 

e var = mm q E(q) . (6.5) 

In accord with the previous discussion let us neglect the 
term in — aq in E(q), which accounts for the gluon ex- 
change (it is discussed in some detail below). The found 
variational binding energy is shown in Fig. [3Jby two thin 
lines. The thin solid line takes into account the total 
energy from (|5.3[) , while thin dashed line neglects the ra- 
diative correction, i.e. discards the last term (V IiL d) in 
this equation. 

Fig. |3] shows that the two approaches to the problem, 
one based on the perturbation theory, another on the 
variational approach, are close. The discrepancy is be- 
low 21% in all the region shown. This discrepancy does 
not contradict the accuracy of the quantum mechanical 
calculation ~ 15%. It also complies with the accuracy of 
the variational approach, which in the problems of this 
type should better than ~ 20%. It is interesting that the 
two approaches also agree on the corrections. Compare 
the radiative correction, which is shown by the discrep- 
ancy of the solid and dashed lines, (the two lines are 
thick for quantum mechanical and thin for variational 
calculation). The observed agreement of two so differ- 
ently formulated approaches is satisfying. 
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FIG. 5: From bottom to top: thick solid, thin solid, dotted, 
dash-dotted and dashed lines: total binding energy e, the ra- 
diative correction (V ra d), the gluon correction — aq, pure kine- 
matic relativistic correction — p 4 /(4m 3 ), and relativistic cor- 
rection (Kei} respectively, all calculated using the variational 
approximation based on Eq. (|5.3[) : mi = 100 GeV, v = 11 
(compare Fig[6l which illustrates the case of mu = 600 GeV). 

Fig. [5] and Fig. [6] illustrate the role played by different 
perturbations in the problem (Fig. [3] is for = 100 
GeV and Fig. [6] is for mu = 600 GeV) . We start our 
discussion from Fig. [5] The solid line there shows the 
total binding energy, which is same as in Fig. [3J and is re- 
peated here to provide the general scale of energies. The 
dashed and dotted lines show the relativistic and radia- 
tive corrections, which are given by the terms (V^ e i) and 
(Vrad) from Eq. (|5.3p . Note that the relativistic correc- 
tion is positive and very large. Remember that this point 
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FIG. 6: From bottom to top: thick solid, thin solid, dotted, 
dash-dotted and dashed lines: total binding energy e, the ra- 
diative correction (V ra d), the gluon correction — aq, pure kine- 
matic relativistic correction —p 4 /(Am 3 ), and relativistic cor- 
rection (Kei) respectively, all calculated using the variational 
approximation based on Eq. (|5.3p ; rrih = 600 GeV, v — 11 
(compare FiglU which illustrates the case of mn = 100 GeV). 



was an obstacle, which was overcome using the approach 
to the perturbation theory formulated in the previous 
subsection. (Note also that here the proper relativis- 
tic correction (V re i) is considered. It should not be con- 
fused with the residual interaction Vr i,res, which arises 
in the perturbation theory, see Eq. (|6.4p . As was men- 
tioned, the latter is smaller and has the negative sign, see 
Fig. [3J) It is interesting that the major contribution to 
the relativistic correction gives the potential energy (V rc i) 
in Eq. (|5.3p .The simple kinematic correction, — p 4 /(4m 3 ), 
which is accounted for by the term (T rc i) in this equation 
proves to be very small, see the dot-dashed line. 

The radiative correction shown in Fig. [5] proves to 
be negative and small, note the blowing up factor ten 
used to make it visible. To emphasize the point we ex- 
tended data up to very large fermion masses m = 2 TeV 
and found that the correction remains very small even in 
this mass interval. The found minor role of the radiative 
corrections in the two-fermion system complies with our 
previous results for the fermion bags (2(| [53] ■ The situ- 
ation for a large Higgs mass is different. The radiative 
corrections are significantly more important there - see 
Fig. 1 



C. Gluon exchange 

The numerical examples presented previously ne- 
glected the processes with the gluon exchange, which 
is acceptable provided the two heavy leptons are con- 
sidered. However, for two heavy quarks the role of the 
gluon exchange needs to be examined. For simplicity 
let us approximate the impact of the gluon exchange by 
the long-range potential V\ r (r) from Eq. (|2.4j) . Then the 
relative strength and impact of the gluon exchange can 
be estimated by comparing the effective coupling con- 



stant a = 4a s /3 in this potential with the Higgs coupling 
= m 2 / (Attv 2 ) . Taking for estimation a s — a s (Mz) ~ 
0.12 one finds that the mass after which the Higgs ex- 
change dominates is mh-dom ~ 349 GeV. Meanwhile, as 
was mentioned, the critical mass that makes the bound 
state due to the Higgs exchange possible is m cr « 515 
GeV. Hence an at this mass substantially exceeds a s , 
a/tC^cr) ~ 2.2(4a s /3), becoming dominant at larger 
masses. For example, at rn = 1100 GeV they differ by 
an order of magnitude, a/,.(1100 GeV) s» 10 (4a s /3). 
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FIG. 7: The gluon exchange in the binding energy of two 
quarks; solid - Higgs exchange only, the energy e with the 
relativistic and radiative corrections (same as the thin solid 
line in Fig. |3J, dot-dashed - additionally, the gluon contribu- 
tion into the energy is included, dotted - naive Coulomb- type 
estimate of the binding energy produced by gluons; calcula- 
tions are based on the variational approach; mi = 100 GeV, 
v = 11. 

One concludes from the mentioned numbers that the 
gluon exchange plays a role of a perturbation, but the one 
which may give a sizable contribution. In order to calcu- 
late it we use the variational approach based on Eq. (|5.3p . 
in which we include now the contribution of the long- 
range potential (Mr) = — ^a s q (omitted in our previous 
discussion). 

The results found for the binding energy are shown 
in Fig. [7J Several engaging features can be observed 
there. Firstly, the discrepancy between the binding en- 
ergy calculated with and without the gluon exchange 
is quite sizable. Secondly, the relative discrepancy be- 
tween these binding energies is obviously falling with the 
fermion mass, which agrees with the fact that the ratio 
&s/ a h is diminishing. However, their absolute discrep- 
ancy remains approximately constant. Thirdly, it is in- 
teresting that the shift of the binding energy produced 
by the gluon exchange greatly exceeds the typical binding 
energy observed when the Higgs exchange is neglected. 
In the latter case one can roughly estimate the spread 
of binding energies due to the gluon exchange alone in 
the Coulomb- type model e g ps e c = — ^ (| a s ) 2 . Fig. [7| 
shows that this energy is much smaller than the shift of 
the binding energy when the gluon exchange is added to 
the Higgs exchange. 
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There is a clear physical reason behind this phe- 
nomenon. The gluon exchange alone keeps two quarks 
at a typical separation that is comparable with the effec- 
tive Coulomb radius. In contrast, the Higgs exchange is 
able to bring the quarks much closer, which makes the 
gluon exchange more effective. 

A similar physical phenomenon takes place in the re- 
gion to < m cr . Here the Higgs exchange alone is not able 
to produce the bound state, but when combined with 
the gluon exchange it increases the binding energy sig- 
nificantly. 



VII. ANALYTICAL SOLUTION AT SMALL 
HIGGS MASS 

For a wide range of fermion masses discussed in the 
data presented in Figs. [3JT] both the fermion mass and 
typical momenta of quarks in the bound state exceed 
the Higgs mass. Consequently it makes sense to investi- 
gate the limit to, — 0. Moreover, the same data shows 
that the radiative correction and the relativistic kine- 
matic correction are small, see Fig. [3 In the simplest 
approximation we can therefore safely neglect these cor- 
rections. Making these simplifications in Eq. (|5.3p we find 
a clear analytical expression for the energy 



E(q) « E a (q) = a h q ■ 

to 



3 



ah q 



(7.1) 



Remember that the second and third terms here originate 
from the Yukawa-type interaction (|2.2p and relativistic 
correction (13.121) correspondingly. To present the argu- 
ment in the most clear form we neglected in Eq. (|7.3|) the 
potential, which stems from the gluon exchange. (The 
way to account it is discussed at the end of this subsec- 
tion). 

The value of the variational parameter q, which pro- 
vides the minimum of E a (q) is 



Qrain = ( ( 1 + 9c# ) l > % - 1 ) . 

The corresponding value of the energy reads 



(l + 9<4) 3/2 -l 



- 1 



(7.2) 



(7.3) 



Remember that a, — m 2 / (Aitv 2 ) is a function of to. 
Hence, Eq. (|7.3l) provides the binding energy of the two- 
fermion state presenting it as a simple analytical function 
of to. Its asymptotic behavior at small values of ah < 1 
is 



64ir 2 v 4 
At large ah > 1 we find 

~ 1 w 



1 



3m 
32n 2 v 4 



2ttv 2 



(7.4) 



(7.5) 



It is instructive to compare these results with the case 
when the relativistic correction, the last third term in 
Eq. (|7.1j) . is neglected. Then clearly the resulting two 
terms constitute the Coulomb-like problem. Their mini- 
mization gives the conventional Coulomb-type energy 



En — — 



64tt 2 u 4 



(7.6) 



Compare now Eqs. (|7.3p and (|7.6p . The first takes into 
account the repulsive relativistic correction, the second 
neglects it. Eq. (|7.4p shows that at small ah there is lit- 
tle difference between the two cases. However, with the 
increase of the mass the distinction becomes prominent, 
compare (|7.5p and (|7.6I) . Fig. |S] illustrates this discrep- 
ancy. The dashed and solid lines there both represent the 
variational energy, in which the relativistic correction is 
included. The difference is that the solid line is taken 
for rrih — ► 0, when analytical Eq. (|7.3p is valid, while the 
dashed one is for to, = 100 GeV, and was calculated nu- 
merically previously, see the data shown by the dashed 
line in Fig. [3] The dotted line in Fig. [5] shows the energy 
from Eq. (|7.6p . which neglects the relativistic correction. 

From Fig. [3] we see again, this time basing the argu- 
ment on clear analytical results, that the relativistic cor- 
rection is important; it substantially reduces the binding 
energy of the two fermions. As a result the system re- 
mains in the nonrelativistic state for much larger masses, 
then one could have anticipated. To illustrate this point 
let us find the value of to, at which the the energy in 
(|7.3p equals the fermion mass. The condition \E a \ = m 
implies a, rs 4.97, which corresponds to the huge mass 
to « 15 TeV. For masses beyond this limit the relativis- 
tic physics takes over, and the developed in the present 
work approach, which is constructed on the nonrelativis- 
tic basis, is not applicable. Remember though the more 
conservative estimate (|3 . 1 8[) for the mass, which limits 
the applicability of our results at to < 0.8 TeV. The just 
presented arguments make it tempting to contemplate 
the possibility that the theory presented remains valid 
for masses beyond the ~ 1 TeV limit, but in order to 
support this claim one should apply more sophisticated 
methods, possibly the Bethe-Salpeter equation, than the 
ones developed in the present work. 

Fig. [5] shows that the dependence of the binding en- 
ergy on rrih is sizable, but still can be considered as a 
perturbations. There is an interesting implication. Re- 
member that previously it was found also that the gluon 
exchange, the radiative correction and the kinematic cor- 
rection — p 4 /(4to 3 ) are all either very small, or at least 
not pronounced. Hence, we can treat them all as per- 
turbations. In order to develop the first order perturba- 
tion theory we can take q mm , which was found neglect- 
ing these perturbations, and substitute it into the com- 
plete expression for the variational energy, which takes 
all these corrections into consideration. 

Remembering that our variational approximation 
matches well the quantum mechanical calculations (< 
20%), we can state that a reliable estimate for the en- 
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FIG. 8: The role of relativistic correction. Solid line - the 
variational energy from (|7.3|) . which presumes mj, = and 
accounts for the relativistic correction; dashed line - the vari- 
ational energy calculated from (|6.5[) presuming rrih = 100 
GeV and taking the relativistic correction into account (but 
without the radiative correction, same data is shown by the 
thin dashed line in Fig. |3j); dotted line - the variational en- 
ergy (|7.6[) . in which rrih = and the relativistic correction is 
neglected. 

ergy e of the two heavy quarks bound via the Higgs bo- 
son exchange provides a simple, transparent analytical 
formula 



greatly the scale of energies typical to the systems bound 
due to the gluon exchange alone. 

The Higgs exchange also strongly increases the bound 
fcrmion density ip 2 (0) and consequently the creation and 
annihilation widths for the bound state. The widths 
calculation will be performed in a separate publication. 
However, there is no fundamental reason, which would 
make impossible the detection of the bound states con- 
sidered at LHC if the heavy fermions do exist. 

Note that we did not calculate the splitting of the spin 
5=1 and 5 = bound states. The relative value of this 
splitting AE/2m should be significantly smaller then the 
leading relativistic corrections since it does not appear in 
the order v 2 /c 2 in the Higgs exchange interaction, which 
is in difference with systems bound via the interactions 
produces by photon and gluon exchanges. 

Note also that if there exist different heavy fermions 
they all contribute approximately equally to the polariza- 
tion operator and vertex in Eqs.((lT]), (|4.3[) . This unusual 
property may be used to measure the number of heavy 
particles via the Higgs-dependent radiative corrections. 
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E(q m in) , 



(7.7) 



in which E(q) was found in (|5 . 3[) while q m - m is given in 
()7.2j) . A more accurate value of g m i n which takes into 
account rrih and the long-range (gluon) potential may be 
obtained using the results of section [TTt 



9b 



((1 + 96 2 ) 1 / 2 -!) 



(7.8) 



where b = (6th + a), see Eq. (|2.15l) . For rrih = a = this 
result coincides with (|7.2[) . 



VIII. CONCLUSION 

We discuss the binding of two heavy fermions due to 
the Higgs boson exchange. Since the effective coupling 
constant increases with the fermion mass, the relativistic 
corrections are important. We find that they are positive 
and significantly reduce the binding energy. However, 
their influence on the critical mass m cr for the fermion 
binding proves to be insignificant. 

The radiative corrections are found to play a minor 
role for all the interval of energies studied, at least up 
to m = 2 TeV. In the case of heavy charged leptons the 
Higgs exchange dramatically increases the binding energy 
when mass m approaches m cr . A very significant increase 
also happens for quarks when m > m cr . The shift of 
energy, which is produced by the gluon exchange in the 
system bound mostly due to the Higgs exchange, exceeds 



Appendix A: Retardation and relativistic corrections 



Taking Eq. (|3.2[) for the relativistic and radiative cor- 
rection we apply the conventional for these problems 
technique replacing lo 2 in the matrix element of the inter- 
action by the two commutators of the scalar interaction 
vertex with the Hamiltonian (compare e. g. (3fjj Section 
83, dH Section 38; one can also compare the quantum 
approach with the classical description in 33] Section 65) 



; pi | A | Pi) (pil/3 2e - iqri |P2) 



(Al) 



; P 'i 



[i?i,fte iqri ] |pi) (p' 2 \[H 2 ,f3 2 e-^) |p 2 ) 



Here |p^) and |p,), i — 1,2 are the Dirac wave func- 
tions describing the scattering states of the two fermions, 
and Hi = a, ■ pi + f3iin are the Hamiltonians for these 
fermions. 

Calculating the commutators in (IA1[) one finds for the 
matrix element, which describes the scattering 



V(pi,p 2 ,q) 



9h 



(A2) 



(7i-(2pi+q)) (72-(2p 2 -q)) 



(q 2 



,2\2 



In the center of mass reference frame p = pi = pi 
it is reduced to V(q) in Eq. ([3.15[l . The potentials in 
12[) and (|3.15p should be sandwiched between the Dirac 
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spinors describing the two fermions. Fulfilling the Foldy- 
Wouthuysen transformation p. 51) we find that the matrix 
elements of the beta-matrices in the nonrelativistic region 
can be replaced by the following expressions 



Pi 



1 - 



1 



Pi 


pi q 


q 2 


2m 2 


2m 2 


8m 2 


pi 


P2 q 


q 2 


2m 2 


2m 2 


8m 2 



+ i 



, (q x pi) • tri 

4m 2 
(q x p 2 ) ■ <t 2 

4m 2 



(A3) 
(A4) 



In (|3.6[) these expressions are presented in the center of 
mass reference frame. Similarly we find that the matrix 
elements of gamma-matrices can be replaced as follows 



71 • (2pi + q) -> - 

72 ■ (2p 2 - q) - 



q ■ (2pi + q) 
2m 

q • (2p 2 - q) 

2m 



(A5) 
(A6) 



Note the dependence of the matrix elements of gamma- 
matrices on q. Another notable feature is an absence of 
the spin dependence. Generically, the spin variable is cer- 
tainly present in the matrix elements of gamma-matrices 
7j — > — (=Fq-H(2p,-±q)- cr 3 ■)/ (2m), the upper and lower 
signs here are for j = 1 and 2 respectively. However, the 



scalar products with the vectors 2pj ± q, which appear in 
(|A5|) and (|A6|) . eradicate this dependence. This complies 
with the scalar nature of the Higgs boson. It also ex- 
plains why the simplified version of the discussion, which 
is adopted in subsection IIII Al and does not appeal to the 
traditional commutator-based technique (jAll) is able to 
match the results of a more robust approach taken here. 
The only available in the problem spin dependence arises 
from the beta-matrices in Eq. (|A3p . (|A4[) and eventually 
leads to the spin-orbit interaction in Eg. (13.141) . 
We substitute Eqs. l|A"3"|) - (|A"6|) into ([A"2l deriving 



^(pi,P2,q) 



~9l 



1 



1 - 



Pi 



■Pi 



2m 2 



4m 2 



(Pi - Pa) • q . (q X Pi) • <r\ - (q x p 2 ) • er 2 
— i 



2m 2 4m 2 
(g-(2pi + q)) (g-(2p 2 -q) 
4m 2 (q 2 + m 2 ) 2 



(A7) 



After that we adopt the center of mass reference frame 
and perform the Fourier transform over q to the coor- 
dinate representation r. Straightforward though lengthy 
calculations lead to Eqs. (pTT2"j) . ([3TT5]) . and (pTH)) , 
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